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Abstract 

We consider a cosmological scenario in which the expansion of the Universe is dominated by 
phantom dark energy and black holes which condense out of the latter component. The mass 
of black holes decreases via Hawking evaporation and by accretion of phantom fluid but new 
black holes arise continuously whence the overall evolution can be rather complex. We study the 
corresponding dynamical system to unravel this evolution and single out scenarios where the big 
rip singularity does not occur. 
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I. INTRODUCTION 



Phantom dark energy fields are characterized by violating the dominant energy condition, 
p+p > 0. Thereby the conservation equation, p+3H(p+p) = 0, has the striking consequence 
that the energy density increases with expansion . In the simplest case of a constant 
ratio w = p/p one has p oc a~ 3( - 1+w \ where 1 + w < 0, while the scale factor obeys a(t) oc 
(t* — t)~ n with n = — 2/[3(l + w)] and t < t*, being t* the "big rip" time (at which the scale 
factor diverges). However, as we shall show below, the big rip may be avoided if black holes 
are produced out of the phantom fluid at a sufficiently high rate. One may think of three 
different mechanisms by which phantom energy goes to produce black holes. 

{%) In phantom dominated universes more and more energy is continuously pumped into 
any arbitrary spatial three-volume of size R. Hence, the latter will increase as a while the 
mass, M, inside it will go up as as a m+3 -with m — 3 |1 +w\. Therefore, the ratio M/R will 
augment with expansion whenever m > —2 (i.e., w < —1/3, dark energy in general). As a 
consequence of the energy being pumped faster than the volume can expand, the latter will 
eventually contain enough energy to become a black hole. This is bound to occur as soon 
as pi? 3 > R/(2G) (i.e., when R > (2G p) -1 / 2 ). 

(ii) As statistical mechanics tells us, equilibrium thermal fluctuations obey < (SE) 2 >= 
ksCvT 2 , where 5E = E— < E > is the energy fluctuation at any given point of the 
system around the mean value < E >, Cy the system heat capacity at constant volume, kp 
Boltzmann's constant, and < ... > denotes ensemble average {3}. As demonstrated in jj], 
unlike normal matter, dark energy gets hotter with expansion according to the law T oc a~ 3w . 
The latter follows from integrating the temperature evolution equation T/T = —3H(dp/ dp) n 
which, on its turn, can be derived from Gibbs' equation T dS = d(p/n) +pd(l/n) and the 
condition that the entropy be a state function, i.e., d 2 S/(dTdn) = d 2 S/(dndT) -see [4] 
for details. Therefore, it is natural to expect that the phantom fluid will be very hot at 
the time when it starts to dominate the expansion. Since the phantom temperature grows 
unbounded so it will the energy fluctuations as well -a straightforward calculation yields 
< (SE) 2 >oc a~ 3( - 1+2w \ Eventually they will be big enough to collapse and fall within their 
Schwarzschild radius. (This simple mechanism was applied to hot thermal radiation well 
before phantom energy was introduced {5]). 

(Hi) As demonstrated by Gross et al. js] thermal radiation can give rise to a copious 
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production of black holes of mass T _1 whenever T is sufficiently high. In this case, black 
holes nucleate because the small perturbations around the Schwarzschild instanton involve a 
negative mode. This produces an imaginary part of the free energy that can be interpreted 
as an instability of the hot radiation to quantum tunnel into black holes. A simpler (but less 
rigorous) treatment can be found in Ref . 7| • One may speculate that black holes might also 
come into existence by quantum tunnelling of hot phantom fluid similarly as hot radiation 
did in the very early Universe. 

Notice that mechanisms (i) and (Hi) differ from the conventional gravitational collapse 
at zero temperature -the latter does not produce black holes when the fluid has w < —1/3, 
see e.g. js]. 

Obviously, one may wonder whether if phantom can produce particles other than black 
holes. In actual fact, there is no reason why dark energy in general (not just phantom) 
should not be coupled to other forms of matter. However, its coupling to baryonic matter is 
highly constrained by measurements of local gravity [9( but not so to dark matter. In any 
case, we will not consider such possibility because it would introduce an additional variable 
in our system of equations below (Eqs. (jSJ)- ©) and would increase greatly the complexity 
of the analysis. 

The number black holes in the first generation (assuming that all those initially formed 
arise simultaneously) will be N { = {a/R)l = C a^ l+w ^ 2 with C = [2Gp a (1+w) }^ 2 -here 
the subscript zero signals the instant at which phantom dark energy starts to overwhelmingly 
dominate all other forms of energy-, and we may, reasonably, expect that they will constitute 
a pressureless fluid. Later on, more black holes may be formed but if much of the phantom 
energy has gone into black holes, then the second generation will not come instantly (as 
there will be less available phantom energy than at t — 

Further, the black holes will accrete phantom energy and lose mass at a rate M = 



— 167rM 2 <fi 2 regardless of the phantom potential, V(<f>) 10(. (Notice that for scalar phantom 
fields p + p = —<f) 2 )- Additional mass will be lost to Hawking radiation 11]. 
The system of equations governing this complex scenario is 



p bh + 3Hp bh = T p x + IQ-nnM 2 (l + w)p x ~ n j^, (!) 
p x + 3(1 + w)Hp x = -Tp x - 16tt nM 2 (1 + w)p x , (2) 
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OL 

Pi + AH Pi = n Jp > ( 3 ) 

3H 2 = p x + p bh + p 1 , (4) 

M = 1QttM 2 (1 + w)p x -j^ (5) 



(subscripts bh, x, and 7 stand for black hole, phantom, and radiation components, respec- 
tively, and we use units in which 8ttG — c — H — 1). Here, T = constant > denotes the 
rate of black hole formation, n = N/ a 3 the number density of black holes, a a positive con- 
stant, and H = a/a the Hubble function. The first term in ([5]) corresponds to the mass loss 
rate of a single black hole via phantom accretion; the second term to spontaneous Hawking 



radiation 12] . For simplicity, we assume that the black holes solely emit relativistic particles 
(i.e., a fluid with equation of state p 7 = p 7 /3). This explains Eq. (j3J). The second term on 
the right of Eq. (j2J) ensures the energy conservation of the overall phantom plus black hole 
fluid in the process of phantom accretion -the mass loss of black holes must go into phantom 
energy. 

At first glance, depending on the values assumed by the different parameters, T, w, and 
a, two very different outcomes seem possible: (i) a big rip singularity, if the phantom energy 
eventually gets the upper hand; and (ii) a quasi-equilibrium situation between the black hole 
and phantom fluids, if none of them comes to dominate the expansion -in this case the big 
rip cannot be guaranteed since the overall equation of state will not be constant. 

There are five equations and, seemingly, six unknowns, namely, (pbh, Px, P-y, H, n and 
M) but, in actual fact only five (pbh, Px, P-y, a , and M) as n can be written as p^h/M. 

To make things easier, we shall further assume that the black holes are massive enough to 
safely neglect Hawking evaporation and that, initially, there was no radiation present. The 
latter assumption is not much unrealistic, the fast expansion redshifts away radiation and 
dust matter very quickly. (The black- hole fluid also has the equation of state of dust but, 
as said above, black holes are continuously created out of phantom). Thus, we can dispense 
with Eq. ([3]) and the last term on the right hand side of Eqs. ([1]) and (jSJ). 

One may wonder whether the black holes may coalesce leading to bigger black holes and 
produce a big amount of relativistic particles. We believe we can safely ignore this possibility 
since the fast expansion renders the chances of black holes encounters highly unlikely. 
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Accordingly, the system of equations reduces to 



p bh + 3Hp bh = Tp x + 16tt (1 + w)Mp bh p x , (6) 

p x + 3(1 + w)Hp x = -Tp x - 16n M(l + w)p bh p x , (7) 

3H 2 = p x + Pbh , (8) 

M = 16kM 2 (1 + w)p x . (9) 

Thus, we may choose the three unknowns, p b h, p x and M -since H is linked to p b h and p x 
by the constraint Eq. (jSJ). 

Let us assume a phantom dominated universe (i.e., no other energy component enters 
the picture initially). Sooner or later a first generation of black holes will arise; then, the 
question arises: "will new black holes condense out of phantom (second generation) before 
the first generation practically disappears eaten by the phantom fluid and consequently, the 
Universe become forever dominated by a mixture of phantom and black holes, or the black 
holes will disappear before they can dispute the phantom the energy dominance of cosmic 
expansion?" 

To ascertain this we shall apply, in the next Section, the general theory of dynamical 
systems \\J$\ to the above set of equations ©-Q. We shall analyze the corresponding 
critical points at the finite region as well as at infinity. As we will see, due to the existence 
of a critical point in the finite region of the phase portrait (connected with the creation rate, 
T, of new black holes), there exist solutions that instead of ending up at the big rip tend 
asymptotically to a Minkowski spacetime. Hence, in some cases, depending on the initial 
conditions, the big rip singularity can be avoided thanks to the formation of black holes out 
of the phantom fluid. 

Before going any further, it is fair to say that, in reality, no one knows for certain if 
phantom fluids have a place in Nature: they may suffer from quantum instabilities jl^j ]. 
although certain phantom models based in low-energy effective string theory may avoid 



them 151 ] . On the other hand, observationally they are slightly more favored than otherwise 
-though, admittedly, this support has dwindled away in the last couple of years. In view 
of this unsettled situation, we believe worthwhile to explore the main possible consequences 
its actual existence may bring about on cosmic evolution. 



II. DYNAMICAL STUDY 



Here we apply the theory of dynamical systems to analyze the above set of differential 
equations (jH])-©- Using the total density px = Pbh + Px, the system can be recast as 



p x = -3(1 + w)y/pr~p x -Tp x - + w)Mp x (p T - p x ) , (10) 

p T = -3(/4 /2 + wp 1 T /2 p x ) , (11) 
M = 16nM 2 (l + w)p x . (12) 

Its critical points follow from setting p x = p? = M = . 
Three possible situations arise, namely: 

1. M = and p x = 0; 

2. M = and p x ^ 0; 

3. M + and p x = 0. 

In virtue of Eq. (flT]) . the first one implies Pt = 0, i.e., the origin. The second one 
corresponds to 

r 

p T = - W p x and y/pr~ = -^7— r . (13) 

3(1 + w) 

Since both densities must be semi-positive definite, we have that w < — 1 which is consistent 
with the assumption of phantom fluid. For fluids satisfying the dominant energy condition 
this critical point does not exist. Finally, the third case also implies the origin as, by virtue 
of Eq. (TIID, p T = 0. 

A. The critical points at the finite region 

In the finite region, the obvious critical point is the origin (M = p x = px = 0), and, for 
the phantom case (w < —1), the point given by 

r 2 r 2 

9w(l + w) 2 9(1 + w) 2 
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After linearizing the system of equations (|T0l) - (|T2l) . we get 



f(l + w)j*= + 16vrM(l + w)p x \ 5p T 
[3(1 + w)y/pf + T + 16vrM(l + w)(p T - 2p x )} 5p x 
16tt(1 + w)(p T - p x )p x 5M, 



5p T 

SM 



^JpT + W- 



'PT 



5px — 3wy/pr dp. 



32vrM(l + w)p x SM + 16vrM 2 (l + w) 5p x . 



(15) 
(16) 
(17) 



1. The critical point at the origin 



In this case, the system reduces to 



Sp x = -T5p x , 8p T = 5M = . 



Here we have made the reasonable assumption that p x /y/pr = as p x , pt 
point is an attractor since its sole eigenvalue is negative. 



(18) 

0. This critical 



2. The finite critical point 

By imposing Eqs. (TH1) on (|T5l) -( fT7j) and linearizing, we get 

r r i6vrr 4 

opx = -7T- d PT + Q1/1 , — dM ' 



°PT = — 0p T + — 0p x , 



1 + W 

SM = 0. 



1 + W 



(19) 

(20) 
(21) 



The roots of the characteristic equation of this system are: 



A, 



r 



2 1 + w 



[l±V-l-2w] 



(22) 



In view that w < — 1, both eigenvalues are real, one positive and the other negative, i.e., a 
saddle point. 
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B. The critical point at infinity 

A full analysis, in the original three-dimensional system, of the critical point at infinity, is 
considerably hard since one has to embed the system in a four-dimensional space of difficult 
visualization. We will consider, instead, the three two-dimensional systems resulting from 
projecting the original one upon three mutually orthogonal and complementary planes. 

1. The critical point at infinity in the plane (p x ,pr) 

By setting M = the system (fTUl) - (|T2l) gets projected onto the plane (p x , pr), resulting 



x = — 3(1 + w)y/yx — Tx = X(x, y) 
y = -3y 3/2 - 3w y/y x = Y(x, y) , 



(23) 
(24) 



where x = p x and y = px- This plane corresponds to the situation that all black holes are 
formed with the same mass which does not vary with time. 

Let us introduce a new, ancillary, coordinate z, to deal with points at infinity, and consider 
the unit sphere in the three-dimensional space (x, y, z). Then, by defining new coordinates, 
u and v, by 



u 



X 



V 

y = -, 



u 2 + v 2 + z 2 = 1 , 



we can write, 



X(x,y) = — 3(1 + w)y/yx — Tx 
1 



?3/2 



3(1 +w)yftu -Tuz 1/2 



?3/2 



P(u,v,z) , 



Y(x,y) = -3y 3 / 2 -3w^ 



x 



y 3/2 



- 3v 3/2 - 3wv 1/2 u 



x3/2 



Q(u,v, z) 



For the two-dimensional system (I23 l) -(j24l ) one follows 



(25) 



(26) 



(27) 



Ydx + Xdy = . 



Moreover, 



dx 



du dz 

^ 2 

z z 2 



and dy 



dv dz 



z z l 



(28) 



(29) 
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Hence 

Adu + Bdv + Cdz = , (30) 

where 

A = -zQ , B = zP, C = uQ-Pv. (31) 
Now, we construct a new three-dimensional system 

it = Bz-Cv, (32) 
v = Cu-Az, (33) 
z = Av-Bu. (34) 

The region at infinity follows by setting z — 0. This implies 

U Q - vP = , u 2 + v 2 = 1 , (35) 

i.e., 

u — , f = 1 , and u = v = —— . (36) 

These two are the critical points at infinity. 

Let us begin by considering the first one, namely, u — 0, v — 1. To this end we perform 
the transformation, 



* = - , »?=-, (37) 



and obtain the following system, 



£u + ft> = Br]v-Cv, (38) 

u = C^u-Attu, (39) 

■fjv + r)v = Av — B£v . (40) 

Then, the two-dimensional system at infinity is, 



i = -C(l + ?)+B V + Arf, 
t) = A(l + rj 2 )-B^-C^rj. 
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(41) 
(42) 



Upon linearizing around the critical point at infinity (£ = r\ — 0, v — 1) we get, 



n 



3r] . 



(43) 
(44) 



Altogether, this critical point at infinity is a saddle point for w > and a repeller for 
w < 0. 



the old u axis comes to coincide with the new axis, v' . Proceeding as before, we obtain 



Clearly, this point is an attractor for w < —1, a saddle for — 1 < w < 0, and a repeller for 
< w. Notice that the big rip singularity corresponds to this point when w < — 1. 

The top panel of Fig. [1] displays the phase portrait. All solutions start from p T = p bh — > oo 
(i.e., the top point on the vertical = 0). Some of them cannot avoid the big rip 

singularity (i.e., the common point to the circle and the straight line x = y). Those solutions 
that end up at the center of the circle (the Minkowski state, p x = pbh = 0) evade the big 
rip. We remark, by passing, that except for the particular case T = 0, no solution can go 
from the origin to the infinity along the straight line p x = Pt- 

For vanishing T (bottom panel of Fig. [T]) the finite critical point collapses to the critical 
point at the origin, which becomes a saddle point. This represents the usual scenario of a 
system composed of pressureless and phantom fluids in which is implicitly assumed that no 
black holes are produced. 

The corresponding phase portraits when the dominant energy condition is satisfied are 
shown in the top (0 > w > — 1) and bottom {w > 0) panels of Fig. [21 



To study the second critical point, u = v = we perform a clockwise rotation so that 




(45) 



(46) 
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FIG. 1: Phase portraits of the system ()23|) - (j24[) when the dark energy is of phantom type. The top 
panel corresponds to the case of a non-vanishing rate of black hole formation. The finite critical 
point (given by Eq. (|14|) ). a saddle, acts a divider: trajectories at its left, and some passing through 
it, avoid the big rip; trajectories to its right, and some passing through it, end up at the big rip 
-see text. The bottom panel corresponds to the case of no black hole production. In this case, all 
trajectories end up at the big rip. 
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y 





2\ x ~- 







FIG. 2: Phase portrait of the system (|23j )- (124j ). with T = 0, when the dark energy obeys the 
dominant energy condition but fails to fulfill the strong energy condition (top panel). For the 
sake of completeness, the bottom panel displays the situation in which the fluid obeys both energy 
conditions. In both cases, because of the absence of phantom fluid, there is neither black hole 
production nor big rip. 

2. The critical point at infinity in the plane (p x ,M) 

The plane pr = does not belong to the physical region since it requires negative energy 
densities. Anyway, for completeness let us analyze this case. After setting px = 0, we get 

x = —Tx+(l + w)yx 2 = X(x,y), (47) 
y = (1 + w) y 2 x = Y(x,y) (48) 

where x = p x and y = 16nM. 
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Performing the transformations 



u = zx , v = zy , u 2 + v 2 + z 2 = 1 , (49) 

it follows that, 

X(x, y) = ^ j - Tz 2 u + (1 + w)u 2 v J = -^P(m, u, 2) , (50) 
Y(x,y) = ^(1 + w)v 2 u = -^Q(u,v, z) . (51) 
As before, the relation 

- Ydx + Xdy = (52) 

implies 

- zQ du + zP dv + (uQ -vP)dz = 0. (53) 

At infinity, z = 0, the relationship uQ — vP = is identically fulfilled. Hence, all points 
satisfying u 2 + v 2 = 1 are singular. This corresponds to a circle centered at the origin. Since 
the latter is an attractor, all the trajectories emanating from the infinity go to the origin to 
end there. So, the infinity is a repeller. 

3. The plane (p T ,M) 

Upon setting p x = the original system reduces to 

p T = -3 py /2 , M = . (54) 

which admits the simple solution, px oc t -2 . Hence, the solution comes from the infinity to 
the origin along the axis p x = 0. 

III. DISCUSSION 

In this paper, we briefly raised the point that black holes may be produced in phantom 
dominated universes by three different mechanisms: (i) via energy accumulation in any 
given spatial three-volume, (ii) gravitational collapse of huge thermal fluctuations, and {in) 
quantum tunnelling of very dense and hot phantom fluid [6]. In this regard, phantom fluids 
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might be viewed as "black holes factories". However, while these processes look rather 
plausible the corresponding calculations are pending. 

Nevertheless, after accepting that black holes may be produced by any of the sketched 
mechanisms, we studied the dynamical system associated to this scenario. Our main finding 
is that because of the existence of a critical saddle point in the finite region of the plane 
(Px, Pt) -given by Eqs. (fl4"j) - the big rip singularity (a generic feature of phantom-dominated 
universes) is no longer unavoidable. This critical point lies in the region of positive densities 
as w < 0. Moreover, since px > p x and Pt/Px = — w, this point is located in the physical 
region (p^r > Px > 0) when w < —1, that is, for phantom equations of state. 

In the said plane (p x ,pt), there is an attractor at the origin, also the finite critical point 
described above (a saddle), and two critical points at infinity: one situated at the axis px (a 
repeller), and another at the straight line px = p x (an attractor). The finite critical point, 
connected with T > 0, is a divider between the solutions that go from p x — > oo to the origin 
and from p x — > oo to the big rip. Those solutions that pass through this critical point, 
depending on the initial conditions, can either go to the origin (i.e., implying that the big 
rip is avoided), or to the critical point at infinity with p x = pt (i.e., big rip) -see top panel 

of Fig. m 

We are now in conditions to fix more precisely the possible scenario described above. 
Let us consider the relations for the finite critical point, Eqs. f|T4l) . Bearing in mind that 
the Universe is nearly spatially flat (Qt = 1), and conceding that this critical point may 
lie not very far from the present expansion era (the Universe began accelerating recently), 
we find T = — 3\/3(l + w)Hq, and £I x q = —w^ 1 . The Wilkinson microwave anisotropy 
probe (WMAP) [3] data are consistent with w ~ —1.1. Hence, for the big rip to be 
attained the present dark energy density parameter must fulfill Q x q > 0.9. Since WMAP 
indicates Q x q — 0.7, our Universe may well avoid the big rip singularity (modulo w is really 
a constant). This also implies T ~ 0.5 H . 

Our system is a three-dimensional one and obviously there are other dimensions. In the 
plane (p x , M) all points at infinity are singular; projecting this onto the Poincare sphere 
we obtain, in that plane, a singular circle around the origin. But, since the origin is an 
attractor, all the points in this circle are repellers. 

As said above, some solutions can evade the big rip. Obviously, the latter become un- 
avoidable if T = (no black hole production), since in this case the critical point at the 
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finite region coincides with the origin. 

There are also other situations in which the big rip can be avoided. For instance, when 



phantom dark energy corresponds to the generalized Chaplygin gas proposed in [17| , or when 
wormholes intervene isj], or when the curvature scalar gets very large and quantum effects 
become dominant [l9|. Nevertheless, to the best of our knowledge, the present scenario was 
never considered in the literature. 

By contrast, as noted by Barrow [20], there are situations in which finite-time future 
singularities can arise even if the fluid filling the Universe obeys p > and p + 3p > 0, i.e., 
under very mild conditions. We do not consider them here. 

Admittedly, it can be argued that in view of the various simplifying assumptions, our 
treatment is not much realistic. In the first place, the rate F is not expected to be a constant, 
it will likely vary with expansion and depend on quantities like w and M. Secondly, we 
have implicitly considered that all black holes are formed simultaneously with the same 
mass -a flat spectrum. It would be more natural to assume the number of black holes 
produced varies with mass and time. Further, as noted earlier, black hole spontaneous 
radiance should be included. Clearly, these features ought to be incorporated in future, 
more realistic, treatments. Nonetheless, we believe this small, first, step may lead the way 
to more ambitious undertakings. 
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